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Abstract. The present work pursues the investigation of the role of spatial 
asymmetry and irreversibility on the dynamical properties of spin systems. We 
consider the ferromagnetic spherical model with asymmetric linear Langevin 
dynamics. Such an asymmetric dynamics is irreversible, i.e., breaks detailed 
balance, because the principle of action and reaction is violated. The fluctuation- 
dissipation theorem therefore no longer holds. The stationary state is however 
still Gibbsian, i.e., the weights of configurations are given by the Boltzmann factor 
corresponding to the ferromagnetic Hamiltonian. The model is exactly solvable in 
any dimension, enabling an analytical evaluation of time-dependent observables. 
We show the existence of two regimes of violation of the fluctuation-dissipation 
theorem in the noncquilibrium stationary state: a regime of weak violation where 
the stationary fluctuation-dissipation ratio is finite but less than unity and varies 
continuously with the asymmetry, and a regime of strong violation where the 
fluctuation-dissipation ratio vanishes asymptotically. This phenomenon was first 
uncovered in the asymmetric kinetic Ising chain. The present results suggest that 
this novel kind of dynamical transition in nonequilibrium stationary states might 
be quite general. We also perform a systematic analysis of several regimes of 
interest, either stationary or transient, in various dimensions and in the different 
phases of the model. 
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1. Introduction 

In the past decades many studies have been devoted to the dynamics of ferromagnetic 
spin systems. If such a system is quenched from a disordered initial state, it relaxes, 
at high temperature, towards its paramagnetic equilibrium state, whereas, in the low- 
temperature ordered phase, equilibrium is never reached and the system undergoes a 
non-stationary coarsening process [1]. 

Both situations arc illustrated by the exactly solvable Glauber-Ising chain [2]. 
At finite temperature, where the relaxation time is finite, the convergence to the 
paramagnetic equilibrium state can be monitored by the exact determination of the 
two-spin correlation functions, either for equal or unequal times. When the system 
has reached equilibrium, its response to a small perturbation is simply related to 
the two-time correlation function, which is the content of the fluctuation-dissipation 
theorem [2] . In the low-temperature regime, where the relaxation time is diverging, the 
system undergoes coarsening, with scaling and universality properties which have been 
characterized in more recent studies [1, 3, 4, 5, 6]. In higher dimensions, since no exact 
result prevails, the present understanding of the properties of kinetic Ising models 
has been achieved by numerical simulations, scaling arguments and field-theoretical 
methods [1]. 

For kinetic Ising models with non conserved (Glauber) dynamics, a line of 
investigation which has been almost completely overlooked until recently is the role of 
spatial asymmetry in the dynamics. This is in contrast with the case of Ising models 
with conserved (Kawasaki) dynamics in the presence of a bias due to an external field, 
which have been thoroughly investigated, mostly numerically [7]. In the context of 
Glauber dynamics, asymmetry means that the flipping spin is not equally influenced 
by all its neighbours. Such an asymmetric dynamics is therefore irreversible, because 
the principle of action and reaction is violated and detailed balance no longer holds. 
In a recent past, two studies have been devoted to Ising models with non conserved 
asymmetric dynamics [8, 9]. Reference [8] gives examples of rate functions for the 
kinetic Ising model in one and two dimensions with totally asymmetric dynamics 
leading to a Gibbsian stationary state (in the sense that the weights of configurations 
are given by the Boltzmann factor corresponding to the ferromagnetic Hamiltonian) . 
Reference [9] investigates, by numerical means, the possible existence of a phase 
transition for the kinetic Ising model, in dimensions two to five, when the usual heat- 
bath rate is modified by truncating the local field, keeping a subset of influential 
spins only. The dynamics thus obtained neither fulfills detailed balance, nor global 
balance (with respect to the ferromagnetic Hamiltonian), and therefore leads to an 
unknown stationary state, except in low dimension, where the problem can be easily 
analyzed [10]. 

The works [8, 9] motivated the systematic study presented in [10]. In particular, 
the examples given by Kiinsch [8] raised the following questions: To what extent are 
these examples unique? Can they be extended to lattices in dimension higher than 
two? As demonstrated in [10], as long as the dynamics is not totally asymmetric, the 
space of parameters defining the rate function of irreversible Gibbsian Ising models in 
one and two dimensions is large. However, imposing total asymmetry of the dynamics 
yields a unique solution (up to a global time scale), for all the examples considered 
(linear chain, square and triangular lattices). The answer to the second question is 
presumably negative: there is no such Gibbsian irreversible dynamics for the cubic 
lattice and one can argue that there are neither totally asymmetric Gibbsian dynamics 
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for lattices of coordination z > 8 [10]. 

Such Gibbsian asymmetric kinetic models, where the flipping spin is not 
equally influenced by all its neighbours, provide benchmarks to study the physical 
consequences of irreversibility, in particular the properties of the resulting nonequili- 
brium stationary state. For instance, for the Ising chain, a one-parameter family of 
asymmetric dynamics can been shown to keep the exact solvability of the symmetric 
Glauber dynamics, allowing a detailed analysis of its dynamical features [11]. One of 
the most striking results is the existence of two regimes of violation of the fluctua- 
tion-dissipation theorem in the nonequilibrium stationary state: when the asymmetry 
parameter is less than a threshold value, there is weak violation of the fluctuation-dis- 
sipation theorem, with a finite asymptotic fluctuation-dissipation ratio, whereas this 
ratio vanishes above the threshold [11]. At low temperature, the asymptotic fluctua- 
tion-dissipation ratio depends on the strength of the asymmetry [3] , and in particular 
vanishes at long times, as soon as the asymmetry is present, while in the absence of 
asymmetry it takes a universal value equal to 1/2 [3]. For the two-dimensional Ising 
model with the asymmetric dynamics introduced in [10] numerical simulations show 
many features in common with the one-dimensional case, in particular the existence 
of two regimes of violation of the fluctuation-dissipation theorem [12]. Irreversibility 
also implies a non-vanishing entropy production rate in the stationary state, which 
can be exactly computed for irreversible Gibbsian Ising models since the stationary 
measure is known [13, 14]. 

In the present work we pursue the investigation of the role of spatial asymmetry 
and irreversibility on the dynamics of spin systems, by considering the ferromagnetic 
spherical model with asymmetric Langevin dynamics. We follow the thread of our 
earlier investigations of the spherical model with symmetric Langevin dynamics [15] 
and of the asymmetric kinetic Ising chain [11]. The outcomes of the present study turn 
out to share many common features with those of the latter reference. In particular 
we show that the phenomenon described above of two regimes of violation of the fluc- 
tuation-dissipation theorem in the stationary state also exists for the spherical model, 
in any dimension, and thereby acquires a status of greater generality. 

The spherical model has been introduced long ago by Berlin and Kac [16], as 
an attempt to simplify the Ising model. It is also known to be equivalent to the 
n — > oo limit of the O(n) Heisenberg model [17]. Its thermodynamic properties can 
be studied exactly in any dimension. The model however possesses non-trivial critical 
properties [16, 18]. The spherical model has another significant advantage in the 
present context: its time-dependent properties under symmetric Langevin dynamics 
can be investigated by analytical means, in all the phases of the model and in any 
dimension [19, 20, 21, 22, 15]. Other aspects of the dynamics of the spherical model 
have been explored, such as the presence of a constant magnetic field, long-range or 
competing interactions, finite-size and surface effects, a temperature-dependent bath, 
and a correlated initial state [23]. The dynamics of spherical spin-glass models has 
also received much attention, especially in the mean-field geometry [24]. 

The spherical model with asymmetric dynamics has also been considered in 
previous studies. One should first mention the study of the dynamics of a mean- 
field spherical spin glass with random asymmetric bonds [25], inspired by neuronal 
networks with asymmetric synaptic couplings. More recently, Langevin dynamics for 
the ferromagnetic spherical model with an asymmetry has been considered in [26] . The 
latter work was however limited to an evaluation of the entropy production rate in 
the stationary state. Instead, the main focus of the present work is on the dynamical 
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properties of the model, at stationarity or out of stationarity, in particular through an 
investigation of the two-time sector of the dynamics, and especially of the fluctuation- 
dissipation ratio. 

The setup of the paper is as follows. The spatially asymmetric linear Langevin 
dynamics is introduced in section 2, and shown to be characterized by a velocity 
vector V. In section 3 we show why the nonequilibrium stationary state reached by 
this dynamics is Gibbsian, in the sense that the weights of configurations are the 
Boltzmann weights corresponding to the ferromagnetic Hamiltonian. Wc also give a 
reminder of some features and properties necessary for the sequel: critical temperature, 
relaxation time, equal-time correlations. In section 4 we derive general expressions for 
dynamical quantities in the two-time sector. As the dynamics is irreversible for any 
non-zero V, the fluctuation-dissipation theorem relating the two-time correlation and 
response functions no longer holds, even in the stationary state. We derive a few 
properties of the nonequilibrium stationary state in section 5, including the dispersion 
of the fluctuation-dissipation ratio in Fourier space and the entropy production rate. 
We then show that there exist two distinct regimes of violation of the fluctuation-dis- 
sipation theorem: a regime of weak violation at small V, where the stationary fluctu- 
ation-dissipation ratio is finite but less than unity, and a regime of strong violation at 
large V, where the fluctuation-dissipation ratio falls off exponentially with time. This 
is demonstrated in one dimension (section 6) and in higher dimension (section 7). The 
next three sections contain a detailed investigation of the transient properties of the 
model, before the nonequilibrium stationary state sets in. We successively consider 
the low-temperature scaling regime of the one-dimensional model (section 8) and the 
vicinity of the critical point in higher dimension, where a continuum description makes 
the analysis simpler, both in the classical regime (D > 4) in section 9 and in the non- 
classical one (2 < D < 4) in section 10. In all these situations, the effect of an 
asymmetry on the dynamical critical properties of the model bears strong similarities 
with the case of the asymmetric kinetic Ising chain. The fluctuation-dissipation ratio 
has a scaling form in two variables, respectively related to the temperature difference 
T — T c and to the asymmetry parameter V. Finally, the asymmetric dynamics of 
the low-temperature ferromagnetic phase is investigated in section 11; it exhibits two 
successive regimes (beta and alpha), in agreement with the phenomenology of glassy 
physics. The quantitative behavior of the two-time correlation function and of the 
corresponding fluctuation-dissipation ratio is shown to be affected by the irreversibility 
of the dynamics. 

2. The model 

The definition of the ferromagnetic spherical model is well known. Consider a lattice 
in arbitrary dimension D 7 chosen to be hypercubic for simplicity. The spins 5 X , sitting 
at the vertices x, are real variables obeying the spherical constraint 



where N is the number of spins in the system, which is supposed to be finite for the 
time being. The Hamiltonian of the model reads 




(2.1) 



X 




(2.2) 



where the sum runs over pairs of neighboring sites. 
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We first remind the rules of the symmetric Langevin dynamics, before addressing 
the case of asymmetric dynamics, which is the focus of the paper. 



2.1. Reversible (symmetric) Langevin dynamics 

Hereafter we assume that the system is homogeneous, i.e., invariant under spatial 
translations. This holds for a finite sample with periodic boundary conditions, and 
for the infinite system. Furthermore, along the lines of earlier studies of dynamical 
properties, we shall consider the mean spherical model, where the constraint (2.1) is 
imposed on average, becoming thus the mean spherical constraint [27]: 

(5 x (t) 2 ) = 1. (2.3) 

Throughout the following, we consistently use the notations of [15]. 

Symmetric Langevin dynamics is defined by the stochastic differential equation 

^=F x -\(t)S x + r lx (t), (2.4) 

where f] x (t) is a Gaussian white noise with covariance 

(T7x(^y(0>=2TS x , y <5(i-t'), (2.5) 

and X(t) is a Lagrange multiplier ensuring the constraint (2.3) at any time t. We 
choose for convenience to parametrize X(t) as 

X(t) = 2D + z(t). (2.6) 

Finally, F x is the force exerted on S x by the other spins, i.e., 

flOY D 

^x = -^=E( 5 *-e«+Sx + eJ, (2.7) 
x a=l 

where the 2D neighbors of x arc y = x ± e a , and e a is a unit vector in direction 
a = 1, . . . , D. The Langevin equation (2.4) thus reads 

d c D 

= ^(Sx-e a + S x+ea - 2S X ) - z(t)S x + Vx (t). (2.8) 

a=l 

This dynamics is reversible and obeys detailed balance with respect to the Hamilton- 
ian (2.2) at any temperature T; it drives the spherical model to the corresponding 
equilibrium state, at least in the paramagnetic phase (T > T c ). 



2.2. Irreversible (asymmetric) Langevin dynamics 

We define the most general asymmetric Langevin dynamics by dissymetrizing the 
roles of the 'left' neighbors (y = x — e a ) and of the 'right' neighbors (y = x + e a ), 
while keeping the linearity of the Langevin equation. We are thus led to deform the 
expression (2.7) of the force F x into 

D 

F x = £((1 + K)Sx- eo + (1 - V a )S x+ea ), (2.9) 

where V a are arbitrary parameters. From now on, the vector V with components V a 
will be interpreted as a velocity. 
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Asymmetric Langevin dynamics is therefore denned by the linear equation 
d<5 D 

a=l 

- z(t)S x + Tu(t). (2.10) 

This dynamics is irreversible, i.e., does not obey detailed balance. It however 
drives the system to a Gibbsian nonequilibrium stationary state, where the weights of 
spin configurations are the Boltzmann weights associated with the Hamiltonian (2.2) 
at temperature T, irrespective of the velocity V. This will be proved in section 3.1. 



3. One-time observables 

3.1. Equal-time correlation function 

From now on, we consider the thermodynamic limit of an infinite system. We assume 
that the initial state is the infinite-temperature equilibrium state with the mean 
spherical constraint (2.3). The spins 5 X (0) at time t = are therefore independent 
and identically distributed Gaussian variables with unit variance. The linearity of 
the Langevin equation (2.10) ensures that the spins S x (t) remain Gaussian at any 
time t. The only non-trivial one-time observable is therefore the equal-time two-point 
correlation function 

C x _ y (t) = <S x (t)S y (t)>, (3.1) 

where brackets denote an average over the infinite-temperature initial state and over 
thermal histories (realizations of the noise). 

The main result of this section is as follows. The equal-time correlation function, 
whose explicit expression is given in (3.13), does not depend on V. It therefore 
coincides with the equal-time correlation function for symmetric Langevin dynamics. 

Hereafter we only summarize the main properties of this correlation function, 
referring the reader to [15] for more details. We have 

C (t) = 1, (3.2) 

reflecting the mean spherical constraint (2.3), and 

C x (0) = <5 Xi0 , (3.3) 

as a consequence of the statistical independence of the spins in the initial state. 

In order to pursue with the explicit evaluation of the correlation function C x (i), 
we have to first solve the Langevin equation (2.10) for an infinite system. This can be 
readily done in Fourier space. Define the Fourier transform by the formulas 



/ F (q) = E /* c ~ iq ' X < /* = / (£p / F (^ eiq ' : 



where 



f d D q 


r d<zi 


f dq D 


1 (2*) D . 




J-* 2tt 



(3.4) 



(3.5) 



is the normalized integral over the first Brillouin zone. Equation (2.10) thus becomes 

= -P*(q) + z (i)]S F (q, t) + /(q, t), (3.6) 
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D 

fi(q) = 2^(1 -cos q a +iV a sin q a ) (3.7) 

a=l 

and 

( V F (c L ,t)r 1 F (< l ',t'))=2T(2ir) D 5 D (c l + 4)S(t-t'). (3.8) 
The solution to (3.6) reads 

^(q^-e-^*-^ 

X ^ F (q,0)+^d S e^) s +^)r, F (q, S )), (3.9) 

with 

Z(t) = [ dsz(s). (3.10) 
Jo 

The equal-time correlation function C F (q, t) in Fourier space is defined by 

(S F (q,t)S F ( q ',t)) = (2^<5 D (q + q')C F (q,i). (3.11) 

By averaging the expression (3.9) over the white noise r? F (q, t), whose variance is given 
by (3.8), and using the initial condition 

C F (q,0) = l (3.12) 

implied by (3.3), we obtain the expression 

p -2w(q)t /ft \ 
C F (q,t) = (l + 2T J dsc 2 ^ s g(s)j , (3.13) 



where 



d 

w(q) = ^(^(q) + "(-<!)) = 2^(1 - co S<Za ) (3.14) 

a=l 

and 

5 (t)=e 2Z W. (3.15) 

The expression (3.14) for w(q) is independent of the velocity V. We have 
thus completed the proof of the result announced above, namely that the equal- 
time correlation function (3.13) is independent of V. As a consequence, one-time 
observables are identical for symmetric (reversible) and for asymmetric (irreversible) 
dynamics. This striking property relies fundamentally on the linearity of the Langevin 
equation. 

3.2. Relaxation time and phase diagram 

We proceed with the calculation of the relaxation time of the dynamics. The Lagrange 
multiplier z(t) 1 and thereby the function g(t) 1 are yet to be determined from the 
constraint (3.2), i.e., 

/(2^F CF(q ' t) = L (3 ' 16) 
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This identity yields a Volterra integral equation [21, 15] for g(t), 

g(t) = f(t) + 2T f dsf(t-s)g(s), (3.17) 
Jo 

with 

f{t) = J (£p e ~ Mq) * = ^ Ul ^)) D (^r D/2 , (3-18) 

where 

I n (z) = J ^Le zcoaq+inq (3.19) 
is the modified Bessel function. In Laplace space, (3.17) yields 

with 

r°° r (]D n i 

f L ( P )= dtf{t)e-*= ^ 9 (3.21) 

Jo J {^r P + 2w(q) 

The latter function has a branch point at p = 0, with a singular part of the form 

/ s » « (87r)- D / 2 r(l - D/2)^" 2 )/ 2 , (3.22) 

as well as a regular part of the form 

/ r L cg (p) = A X - A 2P + A 3 p 2 + • • • , (3.23) 
where the integrals 

A > - /<sf wb? ,3 - 24 » 

arc convergent for D > 2fc. Equations (3.22) and (3.23) jointly determine the small-p 
behavior of / L (p), as a function of dimension D: 

D<2: f L (p) « (87r)- D / 2 r(l - D/2)p-^- D / 2 \ 

2 < £) < 4 : / L (p) (87r)- Z3 / 2 |r(l - D/2)\p^ D ~ 2 ^ 2 , 

D>4: f h (p)^A 1 -A 2 p. (3.25) 

Whenever D = 2, 4, ... is an even integer, the regular and singular parts merge, giving 
rise to logarithmic corrections. 

As is well known, in low enough dimension (D < 2), the spherical model has 
no phase transition. In the present framework, the existence of a phase transition is 
reflected by the divergence of the relaxation time. Such a divergence does not occur 
in low dimension. Indeed f L (p) diverges as p — > 0. As a consequence, at any non-zero 
temperature T, g L (p) has a pole at p = 2p 2 , where T and fj, 2 are related by 

i=2/W) = /(0^. (3.26) 

We have 

9(t) ~e 2 ^, (3.27) 
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r cq = (3.28) 



and so z(t) tends to the constant z(oo) = [i 2 . The system relaxes exponentially fast 
to its noncquilibrium stationary state, with the finite relaxation timej 

1 

which is the same as for symmetric dynamics. The parameter \x defined by (3.26) 
will be interpreted as the inverse correlation length in the continuum limit (see (5.1)). 
The relation (3.28) therefore demonstrates that the dynamical exponent of the model 
has the classical value z = 2, characteristic of diffusion processes, irrespectively of 
dimension. 

The parameter /i is an increasing function of temperature, diverging as 

[i 2 w T-2D + --- (3.29) 

at high temperature (in any dimension D), as can be seen by expanding (3.26), and 
vanishing as T -> 0+ for D < 2, and as T -» T+ for D > 2 (see below). 
On the one-dimensional chain, we have 

fL{p) = ^rrrv (3 - 30) 

so that the dependence of /i 2 on temperature reads 

T = mvV + 4, i.e., n 2 = VT2+4-2, (3.31) 

thus 

A» ~ % (3-32) 

vanishes linearly at low temperature. 

On the two-dimensional square lattice, 

f L ( P ) = k( r^—) , (3.33) 

can be expressed in terms of the complete elliptic integral K. We have therefore 

t = 7t(m 2 + 4) K {wn) A4i k ^' (3 - 34) 

hence: 

H w 4^c~ 27r/T (3.35) 

vanishes exponentially fast at low temperature. 

/n high enough dimension (D > 2), the spherical model has a phase transition at 
a non-zero critical temperature 

Tc=^r=(l^4-X i ■ (3-36) 
2A l \J (2n) D u(tCJ V ; 

Indeed / L (0) = A\ is finite, and so the pole of g h (p) hits the tip of the cut of f h (p) 
at p = at a finite critical temperature T c . The numerical values of T c in various 
dimensions are given in table 1 (section 5). 

X Throughout this work the subscript 'eq' denotes the value of a quantity in the equilibrium state 
generated by symmetric dynamics. 
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Equations (3.26)-(3.28) still hold in the paramagnetic phase (T > T c ). The 
parameter /i is still an increasing function of temperature, diverging as (3.29) at high 
temperature, and vanishing as the critical point is approached from above (T — > T+), 
in a dimension-dependent way, according to 

1/(15-2) 

/ : -I TV \ ~' ■'' ~~ I / --/.. i \ 

2 < D < 4 : n 



|r(i - r>/2)|i? 



L> > 4 



/i 



T — T r 



1/2 



(3.37) 



These estimates obey the scaling law fj, ~ (T — T c )' y , where the critical exponent of 
the correlation length reads 

1 



2 < £> < 4 : i/ = 



D > 4 



D - 2' 
1 



At the critical point (T = T c ), we have 

2<D<4: git) ss (£> - 2)(8^) (I) - 2)/2 sin [(£> - 2)£ 
t— >oo L 2 

In the ferromagnetic phase (T < T c ), we have 

/(t) ^ (STTt)-^/ 2 



9(t) 



M 4 

J "eq 



where the spontaneous magnetization M eq is given by [18] 



Ml 



1 



cq - ^ 

The function g(t) obeys the sum rule 



I 



(3.38) 

t -(4-£>)/2 
7^2 ' 

(3.39) 
(3.40) 
(3.41) 
(3.42) 



4. Two-time observables 

In this section, devoted to two-time quantities, we derive the general expressions (4.3) 
and (4.10) of the two-time correlation and response functions in Fourier space, and 
remind the definition of the so-called fluctuation-dissipation ratio (4.15). 



4-1- Correlation function 

The two-time correlation function is defined as 

C X -y(t,8) = (S X (t)S y (8)), (4.1) 

with < s (waiting time) < t (observation time). 

Its Fourier transform C F (q, t, s) is defined in analogy with (3.11). Using (3.9), 
we obtain 

-Q(q)t-f2(-q)s / rs 
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C F (q, t, s) = C F (q, s) e- n W(*-«) ^j, (4.3) 

where C F (q, s) is given by (3.13). 

In the following, we shall be mostly interested in the local two-time correlation 
function 

C(t, s) = C (t, s) = (S x (t)S x (s)) = J C F (q, t, a), (4.4) 
also referred to as autocorrelation or diagonal correlation function. 

4-2. Response function 

Suppose now that the system is subjected to a small magnetic field H x (t), depending 
on the site x and on time t > in an arbitrary fashion. This amounts to adding to 
the ferromagnetic Hamiltonian (2.2) a time-dependent perturbation of the form 

6H(t) = -Y,H x {t)S x (t). (4.5) 

X 

Within the present framework, it is natural to modify the Langevin equation (2.10) 
by adding the magnetic field H x (t) to the force F x acting on spin S x , so that the 
perturbed equation reads 

D 

= ^((1 + V )Sx-e + (1 - V a )S x+ea - 2S X ) 
a=l 

- z(t)S x + H x (t) + r, x (t). (4.6) 
The solution to this equation reads, in Fourier space, 
^(q^Hc-^)*-^ 



S F (q, 0) + f ds C ^+ Z (s) [H F (q, s) + V F (q, *)] 
Jo 



(4.7) 

It can indeed be checked that the Lagrange multiplier X(t) remains unchanged, to first 
order in the magnetic field. 

Causality and invariancc under spatial translations imply that we have, to first 
order in the magnetic field H x (t), 

(S x (t)) « f ds i?x-y(t, s)H y (s), (4.8) 
Jo y 

where 

S(S x (t)) 



R X -y(t,s) 



5H y (s) 



(4.9) 

is the two-time response function of the model. As a consequence of (4.7), this response 
function reads, in Fourier space 



R F (d,t,s) = e-^-^J 9 ^ ) . (4.10) 
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By (4.3) and (4.10) we obtain the identity 

C F (q,M) = C F (q, S )i? F (q,M) (4.11) 

between the correlation and response functions. Here again, we shall be mostly inte- 
rested in the local two-time response function 

s(s x {t))\ ■ 



R(t,s) = R Q (t,s) = 



SH x (s) 



{ff x (t)=o} 



J^ RF{q ^ s) ' (412) 



also referred to as autoresponsc or diagonal response function. 



4-3. Fluctuation-dissipation ratio 

In an arbitrary stationary state (not necessarily an equilibrium one) , time translational 
invariance implies that two-time quantities only depend on the time difference r = t—s. 
We thus have in particular§ 

C sta t (r) = lim C(s + t,s), # stat (r) = lim R(s + r, s). (4.13) 

s— >oo s— >oo 

Consider now an equilibrium state, i.e., a stationary state for a reversible dyn- 
amics. In such a circumstance, the equilibrium correlation and response functions, 
denoted as C cq (r) and R eq (T), are related to each other by the fluctuation-dissipation 
theorem [28] (see [29] for a general introduction): 



TR cq (r) = - dC ^. (4.14) 



One possible characterization of the distance to equilibrium is the so-called fluc- 
tuation-dissipation ratio [30] (see [31, 32] for reviews): 

ds 

In a nonequilibrium stationary state, the fluctuation-dissipation theorem does not 
hold in general, and so the fluctuation-dissipation ratio remains a non-trivial function 

^stat(r) = lim X(s + r,s). (4.16) 

Likewise, in order to characterize the transient regime before the stationary state 
sets in, we consider the asymptotic fluctuation-dissipation ratio 

X as (s) = lim X(s + t,s). (4.17) 



5. Nonequilibrium stationary state: general properties 

This section as well as the next two ones are devoted to the study of the nonequilibrium 
stationary state attained by the system in its paramagnetic phase, i.e., for T > T c . We 
recall that the critical temperature T c is non-zero for D > 2 only. In this section we 
study general properties of the nonequilibrium stationary state, whereas more specific 
features will be successively discussed in section 6 in the one-dimensional case and in 
section 7 in higher dimension. 

§ Throughout this work the subscript 'stat' denotes the value of a quantity in the nonequilibrium 
stationary state. 
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5.1. Dispersive fluctuation-dissipation ratio 

In the stationary state, the function g(t) grows exponentially as e 2 ^ 2 * (see (3.27)), 
where /i 2 is given by (3.26). The expression (3.13) for the equal-time correlation 
function in Fourier space thus yields the static structure factor 

Equations (4.3) and (4.10) then simplify to 

C F fa T ) - - e -(^(q)+M 2 )r 

iSat(q,T) = e- (n(q)+ " a)T . (5.2) 

These expressions fulfil the stationary form of the identity (4.11), i.e., 

C s F tat (q, r) = C F q (q)i? s F ta t(q, r). (5.3) 

We thus obtain a generalized fluctuation-dissipation relation in Fourier space, of the 
form 

Ti? s F tat (q,T) = _X stat (q) dC£ ^ q ' T) , (5.4) 
where the fluctuation-dissipation ratio 
v , \ w(q) +m 2 

Xstat(q)= 0(q)T^ (5 ' 5) 
(see (3.7), (3.14)) is complex and dispersive, i.e., q-dependent, but independent of t\\. 

At the critical point (T = T c ), all the above results still hold, albeit with \i = 0. 
In the ferromagnetic phase (T < T c ), the static structure factor includes a delta peak 
at q = 0, proportional to the square spontaneous magnetization M 2 q , i.e., 

C c F q (q) - + M 2 q (2^)^^(q). (5.6) 

The expression (3.41) of M 2 q can be recovered by integrating (5.6) over q, using (3.16) 
and (3.36). 



5.2. Entropy production rate 

The evaluation of the entropy production rate in the nonequilibrium stationary state of 
the spherical model with asymmetric Langevin dynamics was the subject of [26]. With 
the notations of the present work, the result for the stationary entropy production 
rate (j> per spin and per unit time reads 

<I>=^(F X (F X -F X )), (5.7) 

where F x and F x are the forces acting on spin for symmetric dynamics (see (2.7)) 
and for asymmetric dynamics (see (2.9)), respectively. Using the symmetries of the 
hypercubic lattice, one can recast the above expression into the following one, 

^((^x-fx) 2 ), (5.8) 

which is more appealing, as it is manifestly positive. 



|| A fluctuation-dissipation relation in momentum space was already considered in [32], 
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The above expression can be further reduced to 

^ = |(i-c 2 , cq )v 2 , 

where 



C, 



2,eq 



T 



A u c\ cos 2qi 



(5.9) 



(5.10) 



(2ir) D w(q) + it 2 

is the static (i.e., equilibrium) correlation function two sites apart along any axis of 
the lattice (e.g. x = 2ei). The stationary entropy production rate only depends on 
temperature and on the square of the velocity vector. 

At the critical point (T = T c , i.e., /x = 0), the expression (5.9) becomes 

^ c -2/ c V 2 , (5.11) 

where 

— cos 2qi 



f d P q 1 



w(q) 



(5.12) 



The result (5.11) also holds for D = 1 and D = 2 at T = T c = 0. Table 1 gives 
the values of the critical temperature T c (see (3.36)) and of the integral I c for various 
values of the spatial dimension D. 



D 


T 

± c 


Ic 


1 





1 


2 





1 - f w 0.363380 


3 


3.956776 


0.209841 


4 


6.454386 


0.146687 


D» 1 


2D 


V(2i?) 



Table 1. Numerical values of the critical temperature T c of the model 
(see (3.36)) and of the integral 7 C entering the expression (5.11) of 
the entropy production rate at T c , for the first few integer values of 
dimension D. Last line: leading asymptotic behavior at large D. 



The entropy production rate remains constant and equal to <f> c in the whole 
ferromagnetic phase (T < T c ). This property is a consequence of (5.6), which implies 
that the difference 1 — C X!C q is proportional to 1 — M^ q = T/T c for all temperatures 
T <T C . The entropy production rate then decreases monotonically from (p c to zero in 
the paramagnetic phase (T > T c ). At high temperature, we have (f> w 2V 2 /T, since 
C2, cq becomes negligible, as it falls off as 1/T 2 . 



6. Nonequilibrium stationary state: one dimension 

Let us now discuss specific features of the nonequilibrium stationary state on the 
one-dimensional chain. The expression (5.1) of the static structure factor reads 



Using the identity 



+oo . , „ 

\p e -iX\x\-iqx = S1Iln ^ (g 2 ) 

cosh (i ~ cos q ' 



Asymmetric Langevin dynamics for the ferromagnetic spherical model 



15 



we obtain the static correlation function 

C x , eri = {S x So)=e-^, (6.3) 
where the inverse equilibrium correlation length 

A = 7— (6-4) 

Scq 

is related to /z and T by 

it = 2sinh^, T = 2sinh/i (6.5) 

(see (3.31)). The expression (5.9) for the entropy production rate simplifies to 

<f> = 2e~i 1 V 2 . (6.6) 

In the stationary state the expressions (4.4) and (4.12) for the correlation and 
response functions become 

2s f do 2(cos<j-iVsing)T 

C stat (r)=T C -(^ / ^ — - 2 

J 2tt 2(1 — cosq) + fi z 

poo 

= T J duc-^+^+^Io (2^(t + u) 2 -V 2 t 2 ) , (6.7) 

= e- (2+ ^ 2)T / (2^1 - V 2 t) , (6.8) 

where Iq is the modified Bessel function (see (3.19)). 

These expressions call for the following comment. The velocity V, introduced 
as a deformation parameter in (2.9), may a priori take arbitrary values. As already 
mentioned, we can restrict ourselves to the case where V is positive. The value V = 1 
turns out to play a special role, demarcating a usual regime (V < 1), where the 
correlation and response functions are positive and decay monotonically to zero, and an 
unusual one (V > 1), where the correlation and response functions exhibit oscillations 
as a function of time. The same phenomenon appears in the asymmetric dynamics of 
the Ising chain [33]. 

6.1. Low-temperature scaling regime 

We first consider the scaling regime where the temperature T rts 2[i w 2/t and the 
velocity V are simultaneously small. In this regime, it is legitimate to treat fi, V 
and q on the same footing, and to consistently use cj(q) rts q 2 and Cl(q) w q 2 + 2iVq. 
This approximation amounts to replacing lattice propagators (given by products of 
Bessel functions) by the much simpler continuum ones (given by Gaussian functions) . 
In the following we shall meet other instances of scaling regimes which lend themselves 
to such a continuum treatment. 

Within this framework, the expressions (6.7), (6.8) for the correlation and 
response functions simplify to 

C st at(r) » \ Y» VT erfc ((/i + V)yft) + c- 2 ^ erfc ((/i -V)y/?)], 

-( M 2 + V 2 )T 



Asymmetric Langevin dynamics for the ferromagnetic spherical model 16 
where 

2 f z 2 

erf z = —j= I dye y , erfcz = l — erf 2 (6.10) 
V 71 " Jo 

are respectively the error function and the complementary error function. The expres- 
sions (6.9) are even functions of V, as it should be, so that it is sufficient to consider 
the case where V > 0. 

The corresponding fluctuation-dissipation ratio X stat (r) is a non-trivial function 
of r, as expected in a noncquilibrium stationary state, which is given by 

» 1 + Vy/¥Fe^ +v '^ 

^stat(T) 

x e- 2 ^ yT erfc ((fi - V)y/r) - e 2 ^ yr erfc {(n + V)y/r) ■ 

(6.11) 

The long-time behavior of this expression depends on the relative values of V and \i. 

• Regime I (V < fi). In this regime, the correlation function falls off as 

2u c -(v 2 +v 2 )t 

C ^ " ^v 2 ~^T> < 6 - 12 > 

while the response function is given by (6.9). The decay rates ac and a R of the 
correlation and of the response are therefore equal: 

a c = a R = [i 2 + V 2 , (6.13) 

and the fluctuation-dissipation ratio has a finite asymptotic value 

u 2 - V 2 

^stat(oo) = 4-^, (6.14) 

which decreases continuously from unity as V increases, and vanishes as V — > [i. 
Regime I is therefore characterized by a weak violation of the fluctuation-dissipa- 
tion theorem. 

• Regime II (V > fi). In this regime, the correlation function falls off as 

C st at(r)«e- 2 ^, (6.15) 

while the response function is still given by (6.9). The decay rates ac and a R 
are now different: 

a c = 2nV <a R = /i 2 + V 2 , (6.16) 

hence the fluctuation-dissipation ratio decays exponentially fast, with a decay 
rate a R — ac = (V — [i) 2 , and so 

X stat (oo) = 0. (6.17) 

Regime II is therefore characterized by a strong violation of the fluctuation- dissi- 
pation theorem. 

• Borderline situation (V = fi). In this case, the correlation function falls off as 

C Bta t(r) w — 2~ , (6.18) 

and so ac — a R — 2fi 2 . The fluctuation-dissipation ratio however decays very 
slowly, as 

1 . 
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The above discussion is very reminiscent of the case of the Ising chain under 
asymmetric dynamics [11]. It is summarized in figure 1, showing a plot of ^" s tat( T ) 
against the ratio V/[X, for several finite values of the dimensionless product ^ 2 t, as 
well as the asymptotic result (6.14). 




Figure 1. Plot of the asymptotic fluctuation-dissipation ratio X a t a t(T) 
pertaining to the nonequilibrium stationary state on the chain in the low- 
temperature regime (see (6.11)), against the ratio V/ n, for fj, 2 r = 1, 2, 5, 10 
and oo (the latter case corresponding to the asymptotic result (6.14)). 



6.2. Arbitrary finite temperature 

We now turn to the analysis of the nonequilibrium stationary state on the chain for 
an arbitrary finite temperature T. 

The response function i? s t a t( - '") falls off exponentially as i? s tat( r ) ~ e~ aRT , with 

a R = 2(l-y/l-V 2 ^ +/A (6.20) 

This asymptotic result can be derived in two ways: either by using the second line 
of (6.8) and the asymptotic growth Iq(z) ~ e z of the modified Bessel function, or 
by evaluating the integral in the first line of (6.8) by the saddle-point method. The 
saddle-point equation reads tan q = —iV. 

The asymptotic analysis of the correlation function C s t a t( T ) is more intricate. 
Consider the integral in the first line of (6.7). Its asymptotic behavior at large r can 
be a priori governed by two special values of q: the saddle point at q"' such that 
tan = — iV, yielding the decay rate 

a {1) = a R (6.21) 

(see (6.20)), and the zero of the denominator at q( u ^ such that cosq^ = 1 + /i 2 /2, 
yielding the decay rate 

aW = Vfi y/fi 2 + 4 = VT. (6.22) 
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Both and turn out to be purely imaginary. Hence the relevant decay rate 
is that pertaining to the special value of q which is closer to the real axis, i.e., whose 
absolute imaginary part is the smaller. We have q^ = q^ (and so — a^) at 
the borderline velocity 

= + 4 = T = tanh £. ( 6 .23) 

The scenario put forward in the low-temperature scaling regime extends to the 
stationary state at any finite temperature. 

• Regime I (V < V c ). In this regime, is closer to the real axis than q^ n \ hence 

Q C = a Ji = a ( ' ) . (6.24) 
The asymptotic fluctuation-dissipation ratio 

^oo)^ 1 - 1 // 1 ^)^ 2 , (6.25) 

i.e., 

i - Ji - v 2 N\ - v 2 

X stat (oo) = V , y c/ v _^L_ ( 6 .26) 

is consistently equal to the value of (5.5) at the saddle point (q = q^). It 
decreases monotonically from unity to zero as V is increased from to V c . Let 
us note that the expression (6.26) has exactly the same form as in the case of the 
Ising chain with asymmetric dynamics [11]. In the latter case, the relationship 
between V c and temperature is given by V c = \J\ — 7 2 , where 7 = tanh2/T. As 
in (6.23), V c increases from to 1, when T varies between and infinity. 

• Regime II (V > V c ). Now q( n ~) is closer to the real axis than q^\ and so we have 

a c = a {II) < a R = a {1) . (6.27) 

The fluctuation-dissipation ratio falls off exponentially fast to zero, with a decay 
rate — a^ n \ 

• Borderline situation (V = V c ). We have then 

2^ 2 A^ 2 +4) T 2 

ac = aR = v^W = = vrm- (6 - 28) 

6.3. Spatial interpretation of the two regimes of fluctuation- dissipation violation 

The existence of two regimes of fluctuation-dissipation violation, i.e., Regime I of 
weak violation, with a finite X sta>t (oo), and Regime II of strong violation, with a 
vanishing X stat (oo), has the following real-space interpretation. Let us start from the 
identity (5.3), which translates in real space into the convolution formula 

+00 

C^statC?") = ^ Crr-^cq-R^statM- (6.29) 
y=-oo 

The static correlation function C x , cq is given by (6.3), while the space-dependent 
response function -R^stat^) reads 

Rx stat(l~) = e~ (2+ ^ 2)T J (i i e i ^+ 2 (cos«-iVsing)r 

= e-( 2 +^ ([±£) ^ I x ( 2v ^^r) , (6.30) 
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where I x is the modified Bessel function (see (3.19)). A very similar expression of 
the stationary response is found for the Ising chain under asymmetric dynamics [11]. 
Consider the expression (6.29) for the autocorrelation function (x = 0), i.e., 

+oo 

C s tat(T) = ^2 C-y tCq Ry, s ta,t(T). (6.31) 

y=-oo 

The first factor falls off exponentially as C-y.eq = c s/l ( see (6.3)), where the decay 
rate fi is the inverse equilibrium correlation length. The second factor, i?y jS tat (t), 
exhibits, in the regime of long times, an exponential profile of the form 

J?0,stat(T" ) 

characterized by the growth rate v or, equivalently, by the asymmetry length £asy m : 
y=- In — - = - . (6.33) 

^ ^ * Sasym 

The two spatial rates fi and i>, or, equivalently, the two lengths £ oq and ( asym , coincide 
for V = tanh/t = V c (see (6.23)). 

In Regime I, i.e., for V < V c , one has v < ft (i.e., £ oq < £ aS y m )- The product 
C-y,cq R y , st a.t( T ) is therefore maximum at the origin (y = 0), and so the sum over 
sites y entering (6.31) is dominated by values of y around the origin. We thus recover 
the result ac — (xr (see (6.24)). 

In Regime II, i.e., for V > V c , one has v > ft (i.e., £ oq > £ aS ym), and so the sum 
entering (6.31) is dominated by distant values of y. A more precise statement is as 
follows. The asymptotic behavior of the modified Bessel function 

I n (z) ~ exp (^Jn 2 + z 2 - nln - + ^ + — ^ , (6.34) 

which holds whenever n and z are simultaneously large, translates into the asymptotic 
fall-off of the response function 

i?»,Btat(r)~e-*W T (6.35) 
in any direction of the x-t plane such that x = wt, where 
tf (to) = 2 + /i 2 - ^w 2 +4(1 - V 2 ) 

w+^w 2 +A(l-V 2 ) 

+ »ln 27TTT0 • (6 ' 36) 

The large-deviation function $f(w) consistently obeys \I/(0) = ap (see (6.20)), as well 
as the symmetry property ty(w) — ty(—w) = —2i>w. It takes its minimum ^ = for 
w = 2V , so that the response function is maximum near the ballistic point x = 2Vt, 
where it falls off as e _Al T . The asymptotic decay rate of the maximum response is 
thus identical to that of the local response for symmetric dynamics (see (6.20)). 
Inserting the estimate (6.35) into (6.31), we are left with the expression 

ac = mm(^(w) — ftw), (6.37) 

w 

where the minimum has to be attained for a positive value of w. We thus obtain 

w = 2(V cosh/t — sinh/t), ac — 2V smhft. (6.38) 

The first expression yields a positive solution (w > 0) for V > V c — tanh/t, whereas 
the second expression is in agreement with (6.5), (6.22). 
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7. Nonequilibrium stationary state: higher dimension 



We now turn to the analysis of specific features of the nonequilibrium stationary state 
in higher dimension. 

The expressions (4.4) and (4.12) for the stationary correlation and response 
functions read explicitly 



Cstat(l") 



-(Q(q)+M> 



J (27r)^ W (q) + M 2 

/■oo p 

T / du C -( 2D +^ T+ ^ J] J o (2V(r + «) 2 -y a V 

"' a=l 



d D q 



-(J2(q)+M 2 )T 



(2tt) d 

e -(W)^/o(2v^^r) 



(7.1) 



(7.2) 



where Iq is the modified Bessel function (see (3.19)). 

These expressions again demonstrate that the usual situation, where the correla- 
tion and response functions are positive and decay monotonically to zero, corresponds 
to the case where all the components of the velocity obey \V a \ < 1. The asymptotic 
behavior of C sta t (t) and i? s t a t (t) can be studied along the lines of the one-dimensional 
case (see section 6). The asymptotic analysis of C sta t(r) again reveals the existence 
of two regimes. 

• Regime I (V small). In this regime, the decay of C sta t(T) and i? s t a t(r), as given 
by the integral representations in the first lines of (7.1) and (7.2), is dominated 
by the same saddle point defined by tangi^ = —\V a . We have therefore 

ac = a R = a^\ (7.3) 

with 

D 



2 ^(i-yrw)+, 2 



(7.4) 



The asymptotic fluctuation-dissipation ratio is given by the value of (5.5) at the 
saddle point q^\ i.e., 



t(oo) 



2E(i-i/y^)+^ 2 



q=l 
D 



(7.5) 



a=l 



This formula is a generalization of the one-dimensional result (6.25). 

Regime II (V large). In this regime, the decay of C s t a t(r) is obtained by constrai- 
ning the saddle point to belong to the manifold with equation w(q) + /J 2 = 0. In 
other words, this decay is dominated by the wavevector q^' so that f2(q( /7 ') is 
extremal while w(q( /7 )) + [i 2 vanishes. We have again 

olc = < a R = a {1 \ (7.6) 
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and the fluctuation-dissipation ratio falls off exponentially fast. 

More explicitly, introducing a Lagrange multiplier £, we are led to extremalize 

the quantity 

Q(q,f) = ft(q)-Mq). (7.7) 

Introducing for convenience the parameter L = 1/(1 — £), this amounts to 
extremalizing 

D 

^(q,L) = 5^(l-cosg +iLV r „sing ). (7.8) 

o=l 

The stationarity equations read tangi^ = — \LV a . We thus obtain 

a^=2LY , V " , (7.9) 

where the parameter L obeys the consistency condition cj(q( 7/ )) + /i 2 = 0, i.e., 
,i 2 = 2yL^-l . (7.10) 

Regime II holds as long as L < 1, so that the imaginary parts of the qi 11 ^ are 
smaller than those of the . 

• Borderline situation. The borderline situation between Regimes I and II 
corresponds to L = 1. We thus have 

M 2 = 2yf^^-lV (7.11) 

The numerator of the asymptotic fluctuation-dissipation ratio (7.5) consistently 
vanishes. At the borderline, we have 

a c = a R = = a W = 2 £ — YL=. (7.12) 

o=l V 1 V a 

This formula is a generalization of the one-dimensional result (6.28). 

The borderline velocities obeying (7.11) run over a closed borderline surface S 
in the D-dimensional space of velocities V. Regime I (weak violation of the fluctua- 
tion-dissipation theorem), with its finite asymptotic fluctuation-dissipation ratio (7.5), 
holds for velocities inside the surface S, while Regime II (strong violation of the fluc- 
tuation-dissipation theorem) holds outside the surface S. The size of the borderline 
surface S increases with /j,, i.e., with the distance to the critical point. 

• For jj, small, i.e., T just above T c , isotropy in velocity space is restored, since V 
is only involved through its norm |V|. In Regime I, (7.4) and (7.5) indeed 
respectively become 

c^^+V 2 (7.13) 

and 

(7-14) 
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The surface S is thus a small sphere of radius \i. In Regime II, (7.9) and (7.10) 
yield 

L = W\ <7 - 15) 

and 

a (7/) =2/x|V|. (7.16) 

• In the opposite limit (fi large, i.e., T large), the surface S approaches the unit 
cube from inside. 

Figure 2 shows a plot of the borderline curve S for the square lattice in the V1-V2 
plane, with equation 

1 1 u 2 

+ ^ =2+^- (7.17) 

yi^T? v / I 3 T f 2 

(see (7.11)), for the values of \i given in table 2, together with the corresponding 
temperatures. 




Figure 2. Plot of the borderline curve S for the square lattice in the V1—V2 
plane, for the values of u listed in table 2. The outermost square is the 
unit square. 





0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.8 


1 


1.2 


1.6 


2 


T 


1.559 


1.888 


2.159 


2.410 


2.654 


2.897 


3.399 


3.936 


4.521 


5.860 


7.454 



Table 2. List of the values of the parameter fi for which the borderline curve 
S for the square lattice in shown in figure 2, together with the corresponding 
values of temperature T (see (3.34)). 
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8. Low-temperature scaling regime in one dimension 

This section and the three following ones are devoted to the transient regime of the 
model, before the nonequilibrium stationary state sets in. The duration of this regime 
is measured by the relaxation time r oq = l/(2fi 2 ) (see (3.28)). The latter is finite 
in the paramagnetic phase (T > T c ), whereas the transient regime lasts forever (the 
system is aging) both at the critical point (T = T c ) and in the whole low-temperature 
ferromagnetic phase (T < T c ). 

In this section we consider the scaling behavior of the one-dimensional model at 
low temperature, in the regime where fi and V are small, whereas times s and t are 
large. The critical regime in higher dimension will be analyzed in sections 9 and 10 
for D > 4 and 2 < D < 4 respectively. Finally, the dynamics in the low-temperature 
ferromagnetic phase (T < T c ) will be the subject of section 11. 



8.1. One-time observables 

Let us start the analysis with the functions f(t) and g(t). In the scaling regime, (3.20) 
and (3.22) (or (3.30)) yield 



/L(P) " Wv 9h{p) " (8 ' 1} 



hence 

fit) « -i= , g(t) « ^= + | e 2 ^ (1 + erf ( P V2t)) . (8.2) 
A quantity of interest in this regime is the reduced dynamical susceptibility 

X(t) = £ C x (t) = C F (0,t) = (^1 + 4/i ^ d Sff ( s )J . (8.3) 

This quantity provides a measure of the range of the correlations, i.e., of the typical 
size of the ordered domains at time t. Using (8.2) and working out the integral in (8.3), 
(or alternatively using (8.7) below) we can recast the dynamical susceptibility into the 
scaling form 

XW~XstatF(x), (8.4) 

where Xstat ~ 2//J, is the static susceptibility, and F{x) is the scaling function 

, n (1 + cr f x) 

F ^ = i , f V ^ 

1 ■ e x (1 + err x) 



XV7T 



with a; = \i\p2± = ^/t/r oq . The scaling law (8.4) interpolates between the coarsening 
regime (t <g; r oq ) and the stationary regime (i ^> r oq ). In the first case, i.e., for id, 
the square-root behavior -F(x) w xy 7 ?!- yields the coarsening law 

X (t) w V8^rt. (8.6) 

In the second case, i.e., for x >• 1, the behavior -F(x) w 1 — 2c _a: /(xt/tt) yields an 
exponential convergence of the dynamical susceptibility to its static value Xstat ■ 
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To close, let us mention that the full structure factor C F (q,t) can be derived 
explicitly in the scaling regime. Starting from (3.13) and using again (8.2), we get 



p -2cft / ft 

CF{q ' t] ^ ~gW V 1 + 4M L Sg[s) 

(j, 2 e 2 » 2t (l+erf((j,\/2i)) + qe~ 2c i 2t (q - i/j, erf (iqV2t)) 



(q 2 + H 2 )g(t) - (8 ' 7) 
This expression interpolates continuously between the Gaussian line-shape 

C F {q,t) « VSnte~ 2q2t (8.8) 
in the coarsening regime and the Lorentzian one 

<*<«> - m 

in the stationary state. 



8.2. Two-time observables 

Let us now analyze the two-time correlation and response functions C(t, s) and R(t, s) 
and the corresponding fluctuation-dissipation ratio X(t, s) in the scaling regime. 
Starting from (3.13), (4.3) and (4.10), and performing the Gaussian integrals over q, 
we obtain 

C(t,s) « -. 1 

-V 2 (t-s) 2 /(t+s) f s e -V 2 (t-s) 2 /(t+s-2u) 



x . h4u / du . g(u) , 

V V^itTs) Wo y/4n(t + s-2u) * V ') 



R(t, s) « - J^l, (8.10) 

where the function g(t) is given in (8.2). Many different sub-cases can be analyzed 
from the general expressions (8.10), according to the values of the dimcnsionless 
combinations [i 2 s, fi 2 t, V 2 s and V 2 t. The most interesting of these regimes are as 
follows. 

Stationary regime. This regime is reached when the waiting time s is much larger 
than r eq {p, 2 s 3> 1). The function g(t) obeys the exponential growth (3.27). The 
corresponding two-time quantities have been studied in section 6.1. 

Coarsening regime. This regime corresponds to the zero-temperature limit of the 
reversible dynamics (fi = V = 0). The resulting scale-invariant coarsening dynamics 
has been investigated earlier [15]. The expressions (8.10) simplify to 

ru ^ ( 4^ 

c(M K(^o 



1 /^ 1/4 



R{t > s) ~7m=T)[v • (8 - n) 
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The corresponding fluctuation-dissipation ratio reads 

and especially (see (4.17)) 

X as (s) ~T(- . (8.13) 

Anomalous aging. This regime corresponds to the zero-temperature limit of the 
irreversible dynamics (/z = 0, while V ^ 0). The expressions (8.10) simplify to 



Us 



1/4 



-V 2 {t-s) 2 /(t+s) 

(t + s) 2 " 
-V 2 (t-s) / f \ !/ 4 

i?(t, S )«— == - . (8.14) 



v^(F- 



The response function ii(£, s) exhibits an exponential decay in r = i — s, characterized 
by the microscopic time scale 1/V 2 . The correlation function C(t, s) however exhibits 
a richer behavior in the late-time regime (V 2 s » 1), with an intermediate anomalous 
aging regime [34] for r <C s, where it assumes the Gaussian form 

C(M) « e- vV2 /< 2s \ (8.15) 

which falls over a characteristic time growing as y/s/V. When r becomes of the order 
of the waiting time s, the correlation function is already exponentially small in the 
variable V 2 s. Its subsequent decay to zero follows an exponential law, characterized 
by the same microscopic time scale 1/V 2 as the response function. The same regime 
can be observed in the case of the Ising chain with asymmetric dynamics [11]. 

Asymptotic regime. This regime corresponds to the limit of a very long time difference 
r = t — s, while the combinations /i 2 s and V 2 s remain finite. It is the non-stationary 
counterpart of the large-r stationary regime investigated in section 6. Here again, 
rather than using the expressions (8.10), it is advantageous to return to the integral 
expressions (4.4) and (4.12). For large r, the latter integrals are dominated by a saddle 
point at q"' = —iV. We thus obtain the following estimates: 



e -(f» 2 +V 2 )r 2 a ( s ) ^ 

C(t, s) « e-"\/^ C F (-iV, s), 



R(t,s)* e e-"\/gW. (8.16) 

The structure factor C F (— iV, s) entering the first of these expressions is given by the 
analytic continuation of (8.7). 

One of the most physically relevant quantities characterizing the transient regime, 
before the stationary state sets in, is the asymptotic fluctuation-dissipation ratio 
^as(s) (see (4.17)). Defining the scaling variables x and y by 



x = ^V2s, y = VV2s, (8.17) 
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the expressions (8.16), together with (8.2) and (8.7), yield the scaling formula 
1 _.o, N (*>v) 



2 ' 



(8.18) 



with 



N(x, y) = ( 1 + 2 (V x 2 ) + 2 ^ 3y2 _ x 2 )e ^ {l + erf x) 



x {x 2 e x \l + ev{x)-ye v \y + xeviy)^ . (8.19) 

In the borderline situation where x — y, i.e., V — V c = fx, the above expression 
simplifies to 

1 „ o , N 0W 



Xas ^ 4x (t¥ + x e x2 (l + erf a;)) 



2 • 



(8.20) 



with 



1 j_ 4 T 2 2 
N (x)= ( V + Ax ° X (1 + erfx) 

^= + 2{x 2 + 1) c* 2 + (2a; 2 + 1) e x * erf x^j . (8.21) 

At the coarsening end, i.e., when x and y are small, the fluctuation-dissipation 
ratio has the series expansion 

- .xy 2 + • • • (8.22) 

The first term agrees with (8.13), whereas only the last one shows a dependence on V. 

At the stationary end, i.e., when x and y are large, the distinction between the 
stationary Regimes I and II shows up as follows. 



• For y < x, i.e., V < jj,, we have 



f_ A y(x + y)(3y 2 - x 2 ) ( ^-^) 



The limit value of this expression is in agreement with the stationary ratio 
^stat(°°) hi Regime I, given by (6.14). In other words, we have 

lim X as (s) = lim X stat (r). (8.24) 

S— ^OO T— ^OO 

This identity testifies that the convergence of the fluctuation-dissipation ratio 
to the value X sta t(oo) is uniform in the s—t plane, irrespective of the direction, 
provided both times s and t are much larger than r oq and 1/V 2 . 
The fluctuation-dissipation ratio X as (s) increases monotonically and reaches 
its limit from below for y/x < l/y/3, i.e., 1/2 < A" as (oo) < 1, whereas it 
exhibits a maximum and reaches its limit from above for l/y/3 < y/x < 1, 
i.e., X as (oo) < 1/2. This feature will also hold in higher dimension, where the 
value X as (oo) — 1/2 will always be singled out. 
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• For y > x, i.e., V > ji, the fluctuation-dissipation ratio falls off exponentially as 

2x 2 (y — x) 



*as(s) 



y(x + y)(3y 2 - x 2 ) 

This corresponds to the stationary Regime II. 

• In the borderline situation where x = y, i.e., V = V c 
pation ratio, given by (8.20), falls off as a power law: 

1 1 
2a; 2 ~ 4/x 2 s' 



(8.25) 

//, the fluctuation-dissi- 
(8.26) 



Figure 3 shows a plot of the asymptotic fluctuation-dissipation ratio X as (s) 
against the scaling variable x = fi\/2s, for several values of the ratio y/x — V/fi. 




Figure 3. Plot of the asymptotic fluctuation-dissipation ratio X as (s) in 
the low-temperature scaling regime in one dimension (see (8.18)) against 
the scaling variable x = ^v2s> for the following values of the ratio V/fi 
(top to bottom): (reversible), 0.2, 0.3, 0.4, 0.5, 1/V3 w 0.577350 (so that 
X sta ,t(oo) = 1/2), 0.7, 0.8, 0.9, 1 (borderline case), 1.2, 1.5 and 2. 



9. Classical critical regime (D > 4) 

In this section we study the transient properties of the model in the vicinity of its 
critical temperature for D > 4, where the critical exponent v = 1/2 assumes its mean- 
field or 'classical' value. We assume that fi and V are small, whereas times s and t 
are large. The dynamical properties of the model obey scaling laws throughout the 
scaling regime thus defined. 
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9.1. One-time observables 

The constants A\ and A 2 denned in (3.24) are finite. As a consequence, (3.20) 
and (3.23) yield in the scaling regime 

hence 

5 (t)«Ge 2 A (9.2) 

The structure factor C F (q, t) can be derived explicitly in the scaling regime. Keeping 
only the integral term in (3.13) and using (9.2), we get 

1 _ c -2( q 2 + Al 2 )* 

c F (q,t)^r c . (9.3) 

The reduced dynamical susceptibility x(t) — C F (0,i) therefore reads 

1 _ „-2 M 2 t 

X(*)«T C - 2 . (9.4) 

M 

This quantity exhibits a linear growth x(t) ~ 2T c i in the coarsening regime, and 
saturates to the static value Xstat = 7c//x 2 in the stationary state. 

9.2. Two-time observables 

In order to study the two-time correlation and response functions, we again start 
from (3.13), (4.3) and (4.10), use (9.2), and perform the Gaussian integrals over q. 
We thus obtain 



C(t,s) w 2T c e-» 2(t - s '> f 
Jo 



e -2^i (s— u)-V 2 (i- S ;r/(t+s-2M) 
Ail 

(4rr(t + a - 2u)) D / 2 



e -(^+v a )(t-«) 



fl(t ' a)W (4^- B ))^ - (9 ' 5) 

The response function does not exhibit any aging property: it coincides with the 
continuum limit of its stationary expression (7.2). Many different sub-cases follow 
from the general expressions (9.5). We again restrict the study to the most interesting 
ones. 

Critical coarsening. This scale-invariant regime corresponds to the reversible dyn- 
amics at the critical point (/z = 0, V = 0) [15]. The expressions (9.5) simplify to 

The corresponding fluctuation-dissipation ratio, 

__. . 1 
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only depends on the ratio t/s. In the limit where this time ratio is very large, the 
fluctuation-dissipation ratio tends to a finite value 

X x = lim lim X(t,s) = I, (9.8) 

which has been emphasized [3, 15] to be a new universal amplitude ratio pertaining 
to nonequilibrium critical dynamics. 



Anomalous aging. This regime corresponds to the irreversible dynamics at the critical 
point (p = 0, while V/0). As in one dimension, the two-point correlation function 
exhibits an intermediate anomalous aging regime for r <C s. The scaling form of 
C(t, s) is however more complex than the Gaussian form (8.15). A scaling analysis of 
the integral entering (9.5), neglecting r with respect to s and integrating over z — s/u, 
indeed yields 



C 21T poo 

C & S ) ~ ( 8 „.)£>/2 s (£>-2)/2 J l dZ C 



V 2 t 2 z/(2s) z (D-4)/2 



T c . „(D-2 



(4tt)^/ 2 (|V|t) b - 2 \ 2 ' 2s 



(9.9) 



where 



T(a,x) = / dyc-Vy"- 1 (9.10) 

J X 

is the complementary incomplete Gamma function. 



Asymptotic regime. This asymptotic regime is the non-stationary counterpart of the 
large-r stationary regime of section 7. The two-time response function (9.5) is already 
in asymptotic form, whereas for large r the integral expression (4.4) for the two- 
time correlation function is dominated by a saddle point at q^) = — iV. Using (9.2) 
and (9.3), we obtain 

_( M 2 +V 2 )r x _ c -2(p 2 -V 2 )s 

^^74^7^ „»-V» • (!m) 
The asymptotic fluctuation-dissipation ratio X as (s) is therefore 

H 2 - V 2 

Xas{s) ~ Ai 2 +V 2 + ( M 2 -3V 2 )c- 2 ^ 2 -v 2 ) S - (9 ' 12) 
Introducing the the scaling variables 

x = liV2s~, y = \V\V2s, (9.13) 
(see (8.17)), the above result reads 

x 2 -y 2 

At the critical end, i.e., when fi 2 s and V 2 s are small, the fluctuation-dissipation 
ratio has the series expansion 

X as (s) « \ (1 + ( M 2 - 3V 2 ) s + ■ ■ •) , (9.15) 

whose first term agrees with (9.8). 

At the stationary end, the distinction between the stationary Regimes I and II 
shows up as follows. 
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• For V 2 < /i 2 , we have 



/i 2 -V 2 ( 1 , 3V 2 -M 2 
^ 2 + V 2 ^ A t 2 + V 2 



* as ( S ) « ( i + ;:, '• ^ v > s ) . o>.i6) 



The limit value is in agreement with the stationary X sta t(oo) (see (7.14)) in 
Regime I. The fluctuation-dissipation ratio increases monotonically and reaches 
its limit from below for V 2 < /i 2 /3, i.e., 1/2 < X sta t(oo) < 1, whereas it 
decreases monotonically and reaches its limit from above for V 2 > /i 2 /3, i.e., 
-^stat(oo) < 1/2- In the special case where V 2 = /i 2 /3, i.e., y — x/y/3, we have 
-^as( s ) = 1/2 for all values of s. 
• For V 2 > /i 2 , the fluctuation-dissipation ratio falls off exponentially as 



V 2 - M 2 
3V 2 - fi 2 



X ^)^^y^^ " )s - (9-17) 



This corresponds to the stationary Regime II. 

• In the borderline situation where V 2 = /i 2 , the fluctuation-dissipation ratio again 
falls off as 

1 1 

2x 2 ~ 4/j, 2 s' 



^-H"^"!^- (9-18) 



Figure 4 shows a plot of the asymptotic fluctuation-dissipation ratio X as (s) 
against x 2 — 2fi 2 s, for several values of the ratio \V\/fi. 




Figure 4. Plot of the asymptotic fluctuation-dissipation ratio X as (s) in 
the classical critical regime (D > 4) (see (9.12)) against x 2 = 2/i 2 s, for the 
following values of the ratio |V| / ju (top to bottom): (reversible), 0.2, 0.3, 
0.4, 0.5, 1/V3 w 0.577350 (so that X as (s) = 1/2 for all s), 0.7, 0.8, 0.9, 1 
(borderline case), 1.2, 1.5 and 2. 
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10. Non-classical critical regime (2 < D < 4) 

In this section we study the transient properties of the model in the vicinity of its 
critical temperature for 2 < D < 4, where the critical exponent v = 1/(D — 2) depends 
continuously on dimension. We again assume that [i and V are small, whereas times s 
and t are large. 

10.1. One-time observables 

Let us, for convenience, parametrize dimension as 

D = 2(l+e) (0<e<l), (10.1) 
so that v = l/(2e). In the scaling regime, (3.20) and (3.23) yield 

n P )^M-Cf, 9\ P )^ pe _ G {2 ^ r (10.2) 



with 



C=(8n)- 1 -e\T(-e)\, G=^. (10.3) 



We have therefore 



g(t)^^(2^) 1 -"h E (2^t), (10.4) 



where the scaling function h e reads^f 



By folding the contour C around the cut (i.e., z — > — p ± iO), we obtain the integral 
representation 



, c e sin7re f°° , e p ^ ,„„„s 
M0 = e*+ / dp . 10.6 

K JO P + P ~ 2 COS 7T£ 

Finally, the series expansion 

trie — 1 j 

TI>1 

relates h £ to the Mittag-LefHer function E £ . 

In the coarsening regime, corresponding to small values of the scaling variable £, 
the scaling function h e diverges as a power law: 

MO^j^r- 1 - (io.8) 

At the stationary end, corresponding to large values of £, it grows exponentially as 

MO*^ (io.9) 

The behavior of the reduced dynamical susceptibility x(t) = C F (0,i) in the 
scaling regime can be estimated by keeping only the integral term in (3.13) and 
using (10.4). We thus obtain 

2T C 



x{t) '" h s (2pH) 



[ dsh £ (2fi 2 s). (10.10) 
Jo 



1 Ci (resp. Co) is a vertical contour in the z-plane such that Rez > 1 (rcsp. < Re 2 < 1). 
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In the coarsening regime (see (10.8)), the susceptibility exhibits the linear growth 
X(t) ~ [2T c t)/e. At the stationary end (see (10.9)), the susceptibility saturates to its 
static value Xstat = T c /(i 2 . 

The following special values of e are of interest, as they correspond to integer 
dimensions. We shall return to them at the end of section 10.2. 

• D = 2 (i.e., e = 0). This case corresponds to the lower critical dimension of the 
model. We have 

ho{0= dx ^=° i+ ^ (10.11) 
Jo r <X> Jo (lnp) 2 +7r 2 

The equality between the two integral expressions is attributed to Ramanujan [35] . 

• D = 3 (i.e., e = 1/2). This case is the only non-trivial one with integer dimension. 
The scaling function reads 

h iMO = ^ + ^Hl + cv^)- (10-12) 

• D = 4 (i.e., e = 1). This case corresponds to the upper critical dimension of the 
model. The scaling function becomes a simple exponential: 

fci(0 = e* (10.13) 

(see the comment at the very end of section 10.2). 

10.2. Two-time observables 

Starting from (3.13), (4.3) and (4.10) and performing the Gaussian integrals over q, 
we are left with the following expressions for the two-time correlation and response 
functions in the scaling regime: 

OT f s p -V 2 (t-s) 2 /(i+s-2«) 

C(t,s) sa -= c = / du- — ; h £ (2ii 2 u), 

In what follows we restrict the study to the most interesting of the different sub-cases 
which can be investigated from these general expressions, according to the values of 
the dimensionless combinations fi 2 s, [i 2 t, V 2 s and V 2 t. 

Critical coarsening. This scale-invariant regime corresponds to the reversible dyn- 
amics at the critical point (fj, = 0, V = 0) [15]. Owing to (10.8), the expressions (10.14) 
simplify to 



C(t,s) 



4T r t ZO- D/4 



(D - 2)(4n) D / 2 (t - s )(D-2)/2 t + 



s \ s 



1 ft^- D/4 



R ^~ {M t-s))^ [s) ' (10 ' 15) 

The corresponding fluctuation-dissipation ratio reads 

X(t,s)Ki ? . (10.16) 

2 ft-s^ 2 



1 + D - 2 \t + s 
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The asymptotic fluctuation-dissipation ratio assumes the universal value [15] 

*co = T ^ L - (10.17) 

Anomalous aging. This regime corresponds to the irreversible dynamics at the critical 
point (p, = 0, while V ^ 0). The two-point correlation function again exhibits an 
intermediate anomalous aging regime for r <C s. A scaling analysis of the integral 
entering (10.14), neglecting r with respect to s and integrating over the dimensionless 
variable z — s/u, yields 

C{ ' j (47r)^/2(|V|r)^2 e • 

Asymptotic regime. The asymptotic regime is the non-stationary counterpart of 
the large-r stationary regime analyzed in section 7. For large r, the integral 
expressions (4.4) and (4.12) are dominated by a saddle point at q^' = — iV. We 
thus obtain the following estimates: 

e -(M 2 +V 2 )r 2 

M (4^ T) ^/2 s 7MV^c F (-iV, s ), 

e -(M 2 +V 2 )r 2 

s ) « (47rT)g / 2 S vW^), (10-19) 

CF (- iV ' s ) « [ S due-^h E (2^u). (10.20) 

h e {2n 2 s) J Q 

The asymptotic fluctuation-dissipation ratio X as (s) thus reads 

1 2(3V 2 M2^ S )-^(2^ S )) 
A as (s) ~ + M 2 M 2 *) 2 

x f due 2V ^ s -^h £ (2^ 2 u). (10.21) 
Jo 



with 



Setting 



s V 2 u 2 

e = * 2 = 2 M 2 S =-, -=^ = ^ ( 10 - 22 ) 

(see (9.13)), we obtain the scaling form 

1 »2 + * h f-?<® f^e^hM. (10.23) 
A as (s) M£) Jo 

At the critical end, i.e., when x and y are small, the fluctuation-dissipation ratio 
has the expansion 

asl ' l + e \ + 2T(2e) + (1 + e) 2 V + )' [ > 

The first term coincides with (10.17), whereas the last one exhibits the leading 
dependence of the result on V. For e = 1, i.e., D = 4, the full expansion agrees 
with (9.15). 

At the stationary end, the distinction between the stationary Regimes I and II 
again shows up as follows. 
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• For V 2 < ^i 2 , the fluctuation-dissipation ratio converges to its stationary value 
X stat (oo) (see (7.14)). 

• For V 2 > /i 2 , it falls off exponentially as 

X as (s) ~ c- 2 ^ 2 -^. (10.25) 

• In the borderline situation where V 2 — n 2 , it again falls off as 

1 

4^i 2 s' 

To close, we analyze in some more detail the asymptotic regime in the following 
cases, where the dimension is an integer. 



« 7-2". (10.26) 



The two-dimensional case (e = 0). This marginal situation corresponds to the lower 
critical dimension of the model. We thus have T c = 0, while n vanishes exponentially 
fast at low temperature, according to (3.35). The critical end is affected by the 
presence of logarithmic corrections to scaling. Let us introduce for convenience the 
logarithmic variable 

A = ln* =ln^- =ln^, (10.27) 
£ 2^ 2 s s y 1 

which is large and positive in the coarsening regime (s <C T cq ). In this regime, the 
rightmost expression for /io(£) m (10.11) is dominated by the integral. Setting u = 
we obtain the estimate 



^"J, dM (A + ln M) 2 + 7 r 2 ' ^ 
which admits the asymptotic expansion 

1 2C 3C 2 — — 

£ MO « ^ + + + ■■ ( 10 - 29 ) 

as A — » oo (i.e., £ — > 0), where C w 0.577215 is Euler's constant. We are thus left 
after some algebra with the logarithmic expansion 



1 , C , , r 2 , 2(4C-3) 



*as(*)~ A + A2+---- V 2 ^X + A 2 +---J- (10 ' 30) 

The asymptotic fluctuation-dissipation ratio thus starts increasing as an inverse 
logarithm of time in the coarsening regime in dimension two, whereas its behavior 
at the stationary end follows the generic pattern, with its two regimes. 

The three-dimensional case (e — 1/2). More explicit results can be derived in 
this case, which is the only non-trivial one with integer dimension. Inserting the 
expression (10.12) of the scaling function hi/ 2 (Q into (10.23), we obtain the following 
scaling form for the fluctuation-dissipation ratio: 

*V 2 + 7^ < 10 ' 31 > 

Aas[S> 2(x 2 -y 2 ) +xc- 2 (l + crfa;)J 

with 

N(x,y)= ( 1 + 2(3 f ~ ^ + 2x{3y 2 - x 2 ) <?* (1 + erf x) 



x [xe x \l + er:ix) - e y2 (x + y erf yf) , (10.32) 
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where the scaling variables x and y have been defined in (9.13). In the borderline 
situation where x = y, the above expression simplifies to 



1 



N (x) 



with 



*as(s) 



N (x) 



Ax f-J= +xc x2 (l + erf x) 



(10.33) 



1 + Ax 2 



+ Ax 6 e x (1 + erfx) 



2x 



2x A e x +(2x 2 + l)e x erfx 



(10.34) 



The above expressions bear a strikingly close resemblance with the results (8.18), (8.20) 
describing the low-temperature scaling regime in one dimension. 

At the critical end, i.e., when x and y are small, the fluctuation-dissipation ratio 
has the series expansion 



^as(s) 



64 - 15tt 



20 



y 2 + ' 



(10.35) 



6 108 27 ' 

The first term agrees with (10.17), i.e., Xao = 1/3 for D = 3. 

Figure 5 shows a plot of the asymptotic fluctuation-dissipation ratio X as (s) 
against the scaling variable x — fi\/2s, for several values of the ratio y/x — |V|//x. 




Figure 5. Plot of the asymptotic fluctuation-dissipation ratio X as (s) in 
the critical regime in dimension three (see (10.31)) against the scaling 
variable x — ^v2s, for the following values of the ratio y/x = |V|//i 
(top to bottom): (reversible), 0.2, 0.3, 0.4, 0.5, 1/V3 w 0.577350 (so that 
X s tat(oo) = 1/2), 0.7, 0.8, 0.9, f (borderline case), 1.2, 1.5 and 2. 
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The four- dimensional case (s = 1). This situation corresponds to the upper critical 
dimension of the model, where the exponent v changes from non-classical to classical 
(see (3.38)), so that logarithmic corrections to scaling can be expected a priori. 
The scaling function is, however, a simple exponential (see (10.13)). As a 

consequence, to leading order in the critical regime, the quantities which can be 
expressed in terms of are given for D = 4 by their expressions in the classical 

regime (D > 4), derived in section 9. In other words, logarithmic corrections to 
scaling are subleading for these quantities. This phenomenon was already put forward 
by Ebbinghaus et al [36]. In particular, the asymptotic fluctuation-dissipation ratio 
is given by (9.12), up to correction terms which are of the order of l/ln(l//u 2 ), i.e., 
1/| ln(T - T c )|, in relative value. 



11. Ferromagnetic phase 

This last section is devoted to the two-time observables of the model in its low- 
temperature ferromagnetic phase (T < T c ). This ordered phase only exists in high 
enough dimension (D > 2). It is characterized by a non-zero value of the spontaneous 
magnetization M eq , given by (3.41). 

The dynamical properties of the model in the late-time coarsening regime are 
ruled by the power-law decay (3.40) of the function g(t). In particular, the behavior 
of the reduced dynamical susceptibility %(i) = C F (0,t) can be estimated by keeping 
only the integral term in (3.13), and using (3.40) and the sum rule (3.42). We thus 
obtain 

X (t) « M 2 c ^t) D '\ (11.1) 

This result has a simple interpretation: the typical linear size of a magnetized domain 
grows proportionally to the diffusive scale L(t) <~ \/t [1], so that its volume, measured 
by x(t), grows as L(t) D ~ t D / 2 . 

The commonly accepted phenomenology of glassy dynamics (see [37] for reviews) 
tells us that two-time quantities behave differently in the following two regimes, which 
will be studied successively hereafter: 

• a beta (stationary) regime (t = t—s <C s), where the two-time correlation function 
C(t, s) only depends on r, and decays from unity to its plateau value M c 2 q . 

• an alpha (aging) regime (r <~ s), where the two-time correlation function further 
decays from its plateau value M 2 to zero. 



11.1. Beta (stationary) regime 

In this first regime, i.e., for r <C s, it is legitimate to simplify (4.3) and (4.10) to 
Cf (a t) w T p-^i) 7 

LJ(q) 

< ta (q,r)«e^(^. (11.2) 

These expressions coincide with the [i — » limit of (5.2), which was shown to also 
describe the nonequilibrium stationary state at the critical point. 
The correlation function therefore behaves as 

Cbeta(T) w M c 2 q + (1 - M c 2 q )<&(r), (11.3) 
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where the function 

f d D a c~° (q)T 

Hr)=T c / q D (11.4) 
J (2-k) u w(q) 

decreases from $(0) = 1 to $(oo) = 0, and thus describes the relaxation of correlations 
throughout the beta regime. 

Hereafter we focus our attention onto the scaling regime where V is small. The 
response function assumes the simple expression 

-V 2 r 

^beta(T) « (4 ^ t)I3/2 , (H.5) 

and therefore falls off exponentially to zero. 

More interestingly, the tail of the function $(r), describing how the correlation 
function tends to its plateau value M 2 q , assumes the scaling form 

*fr> * (4^/vW 2 )/ 2 F (1L6) 
where the scaling variable is 

z = V 2 t, (11.7) 



and with 



where 



F( Z ) = Z -^-^ 1 (^,z), (11.8) 



7(a,x) = [ 
Jo 



Aye- y y a - 1 (11.9) 



is the incomplete Gamma function. The scaling function F(z) starts from the finite 
value F(0) = 2/{D - 2), and decays as F(z) w T((D - 2)/2)z-^ D - 2 ~>/ 2 at large z. As 
a consequence, we have 

2T 

$W ^ (^-2)(4^/M^)/ 2 (1L10) 
for reversible dynamics (V = 0), and 

, . T C T((D -2)12) 

as soon as V ^ 0. 

The correlation function therefore always exhibits a slow power-law convergence 
toward its plateau value. The exponent of this power law however depends on whether 
the dynamics is reversible or not. In the first case (see (11.10)), it equals (D — 2)/2; 
in the second case (see (11.11)), its value is doubled to D — 2. 

Finally, in the same scaling regime where V is small and r is large, the fluctua- 
tion-dissipation ratio also assumes a scaling form, namely 

^beta(r) « X(z), (11.12) 

with 

2e" 2 
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The fluctuation-dissipation ratio is identically equal to unity for reversible dynamics 
(V = 0) throughout the beta regime. In the irreversible case, it starts decreasing as 

4z 

*(*) = 1 - — + .-., (11.14) 

i.e., Xbota(T) ~ 1 — 4V 2 t/.D, whereas it falls off roughly exponentially at large z. 

To close, we give explicit expressions for the scaling functions F(z) and X(z) in 
dimension three and four: 

= 4: F(,) = l=p, 4) ^,,,.',,., (11-15) 
Figure 6 shows a plot of the scaling function X(z) for D = 3 and D = 4. 




Figure 6. Plot of the scaling function A^z) describing the fluctuation- 
dissipation ratio in the beta regime in the ferromagnetic phase, against 
z = V 2 r, for D = 3 and D = 4. 



Alpha (aging) regime 

Let us again consider, for the sake of simplicity, the scaling regime where V is 
small. Throughout the alpha (aging) regime, where both times s and t are large and 
comparable, (3.13), (4.3) and (4.10), together with (3.40) and the sum rule (3.42), 
lead to the estimates 

CaWM) - Me 2 q (j^) D/4 ^ 2{t - S)2/it+S) > 

tfalpha^, S) « -jjj^ - , (11.16) 
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which bear a strong resemblance with the formulas (8.14) pertaining to the zero- 
tcmpcrature dynamics of the one-dimensional model. 

The corresponding fluctuation-dissipation ratio reads 

at ft a- S \( D + 2 )/ 2 

alphal ' ' D(87r)D/m^ \t-s) 

D(t + S) „-2V 2 s(t-s)/(t+s) 



" D(t + s)+4V 2 s(3t + s) ' 1 j 

The first line is the known result for reversible dynamics [15, 22], whereas the second 
one describes the effect of a (small) velocity V in the aging regime. 



12. Discussion 



We have introduced a spatial asymmetry into the linear Langevin dynamics for the 
ferromagnetic spherical model on the hypercubic lattice. The asymmetry is measured 
by an arbitrary velocity vector V. The resulting dynamics is irreversible. It therefore 
breaks detailed balance and its numerous consequences, including the fluctuation- 
dissipation theorem. The corresponding noncquilibrium stationary state is however 
still Gibbsian, with the weights of configurations being dictated by the ferromagnetic 
Hamiltonian. 

The model remains exactly solvable in any dimension, allowing thus an analytical 
evaluation of time-dependent observables. The main emphasis has been put on two- 
time quantities, and especially on the fluctuation-dissipation ratio. We have performed 
a systematic investigation of several regimes of interest, either stationary or transient, 
for arbitrary dimensions and in the different phases of the model. One of the most 
noticeable outcomes of this study is the existence of two regimes of violation of the fluc- 
tuation-dissipation theorem in the nonequilibrium stationary state: a regime of weak 
violation at small V, where the stationary value of the fluctuation-dissipation ratio 
is finite but less than unity, and varies continuously with V, and a regime of strong 
violation at large V, where the fluctuation-dissipation ratio vanishes asymptotically. 
This phenomenon had been uncovered for the first time in the kinetic Ising chain 
under asymmetric dynamics [11]. The present study suggests that this novel kind of 
dynamical transition in nonequilibrium stationary states might be quite general. We 
have also characterized the rounding of the above transition in the transient regime. 

As mentioned in the introduction, numerical simulations of the two-dimensional 
Ising model with the asymmetric dynamics introduced in [10] are currently in 
progress [12]. Finally, it would be worth pursuing the exploration of irreversible 
dynamics and nonequilibrium stationary states driven by spatial asymmetries in other 
directions, such as conserved dynamics, or with other approaches, such as field- 
theoretical methods. 
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